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When fluid is confined between two molecularly smooth surfaces to a few molecular diameters, it 
shows a large enhancement of its viscosity. From experiments it seems clear that the fluid is squeezed 
out layer by layer. A simple solution of the Stokes equation for quasi-two-dimensional confined flow, 
with the assmption of layer-by-layer flow is found. The results presented here correct those in Phys. 
Rev. B, 50, 5590 (1994), and show that both the kinematic viscosity of the confined fluid and 
the coefficient of surface drag can be obtained from the time dependence of the area squeezed out. 
Fitting our solution to the available experimental data gives the value of viscosity which is ~7 orders 
of magnitude higher than that in the bulk. 
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PACS numbers: 68.35.-p, 68.45.-v, 68.15.+C 



The observed giant enhancement of viscosity dealing 
liquids under confinement has led to many experiments 
P* BBB' computer simulations 0, H, > an d some the- 
ory IS IE Oil • When confined to a spacing of few mono- 
layers between two solid surfaces, fluids start behaving 
very differently than in the bulk, exhibiting a dramatic 
increase in viscosity of six to seven orders of magnitude, 
and perhaps an yield stress P, [llj . 

The literature is conflicting on the question whether 
the observed slowing down of the dynamics in the con- 
fined fluid stems from a phase transition 11]. For fluid 
particles of regular shape, some interpretations show the 
possibility of surface induced crystallization [l^, others 
of a glass transition [llj] , both at temperatures well above 
the respective bulk transition temperature. It is some- 
what unsatisfactory, if not strange, that the same basic 
experiment is interpreted in both terms. One important 
conclusion can be drawn from these that our theoretical 
understanding of the phenomenon is not really complete. 

In computer simulations, using surfaces with periodic- 
ity commensurate with the bulk crystalline phase of the 
confined fluid, the fluid showed an ability to sustain an 
applied static stress [12|- But in the simulation of Gao 
et al 14] with incommensurate surfaces, this does not 
happen. This shows the importance of commensurabil- 
ity. But the experimental systems studied so far involve 
confined fluids with a variety of molecular shapes and 
architecture, and are generally incommensurate with the 
surfaces. Therefore commensurability is probably not rel- 
evant |l5| . Actually confinement is a means to introduce 
geometrical constraints [l6| and this slowing down of the 
dynamics should have a more general origin rather than 
a mere commensurability. 

The high increase in the viscosity is perhaps a re- 
sult of mode-coupling [l|| [lS| modified by confinement 
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[r7| . Such a theory is beyond the scope of this work. 
It is however interesting that a simple phenomenological 
treatment in the form of layer-by-layer viscous flow gives 
physically relevant results. 

This model introduced here was first considered by 
Persson and Tosattijg]. But their solution does not say 
anything about the kinematic viscosity, since their cal- 
culation seemed to show that the spreading of the mono- 
layers is determined by the wall-fluid friction, which they 
have introduced as a separate parameter, and not by the 
kinematic viscosity of the fluid. In fact the spreading 
does depend on the viscosity, as shown below. The Pers- 
son and Tosatti model is thus more interesting than its 
authors showed it to be. I redo the Persson and Tosatti 
calculation below with a result different from theirs, and 
I point out where there was an error in their calculation. 
The result is important because it shows that a measure- 
ment of the time evolution of the area of the squeezed 
monolayer gives an estimate of the viscosity. The value 
of viscosity extracted by fitting the time-dependence of 
the area to the form obtained here is consistent with those 
measured in independent experiments [4j . The value of r\ 
(the coefficient of drag by the solid surface on the liquid 
film) is obtained to be of the same order of magnitude as 
that found by Persson and Tosatti 

Let us take the geometry to be circular of radius Tq. 
The plates, confining the fluid are taken to be smooth 
and devoid of any imperfections. Let Pq be the normal 
force applied to squeeze out the fluid. Let us assume that 
when the fluid is confined to a few molecular thickness, 
monomolecular layer of liquid are squeezed out one at a 
time as a result of the applied normal force as shown in 
Fig.l. This crucial assumption, the main phenomenolog- 
ical input in this approach, is supported by the experi- 
ment of Becker and Mugele [2(| • As in the experiment of 
jlj, the upper limit of change of liquid density p due to 
shear is 4%. So we can take the fluid to be incompressible 
in our model. 

Suppose one monolayer is squeezed out, and let v(r, t) 
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FIG. 1: The suggested flow geometry of Persson and Tosatti. 
The squeezing flow results in the nth layer becoming an annu- 
lus with a central hole of radius r, which spreads and pushes 
the rest of that layer outward, while the (n — 1) layers below 
are unaffected. 



be the in-plane velocity field of this monolayer as a func- 
tion of the two-dimensional in-plane coordinate r and 
time t. Then incompressibility and the Navier-Stokes 
equations become: 



V-v(r,t) = 



<9v „ 1 

— +v.Vv = 

dt mn a 



Vp + v V v — r\w 



(1) 



(2) 



where, p is the two-dimensional pressure, v is the kine- 
matic viscosity, rj is the coefficient of drag force by the 
sheets on the fluid, mn a is the two dimensional mass 
density. 

The flow is very slow, so we ignore inertia in Eq. J5J. 
The suggested flow geometry of Persson and Tosatti is 
that the squeezing flow results in the nth layer becoming 
an annulous with a central hole of radius r, which spreads 
and pushes the rest of that layer outward, while the (n — 
1) layer below are unaffected. Then 



v = v(r, i)f 

p = p(r,t). 
Using Eq. (3) in Eq. (1) we get, 

-|-M = o 

r or 



or, 



v dv 

~ + 7T = °- 
r or 



(3) 



(4) 



(5) 



(6) 



Equation (2) becomes 



1 dp , 1 d dv 



dv 

dt mn a dr r dr dr' 

where symmetry implies only the radial component is 
non-zero. 

At this point, Q makes a mistake by making the vis- 
cosity term, i.e. the second term in right hand side of 



Eq. (7), to be equal to zero using Eq. (6). But this term 
should actually be equal to v^. Using this value Eq. (7) 
becomes 



dv 

lit 

From Eq. (6), we get 



1 dp v 

1 v— - rjv 



mn a dr 



r = v 



B(t) 



dv _ B 
~di~~' 



(8) 



(9) 



(10) 



where B(t) is an integration constant, which depends 
only on t. Now using the above two results in Eq. (8) we 
get, 



dp 

dr 



(B + nB) B 

-mn a \-nin a v— (11) 



In solving this equation, let us take the boundary condi- 
tion as p=po at r=ro . Then, 

r mn a vB 1 1 
P(r) -po = -mn a (B + r]B )ln(— ) (-^ - — ) 

ro I r A rg 

(12) 
Now using Eq. (9) and Eq. (10), the above equation can 
be written as 



p(r)-Po = - 



mn a rj d 
4 dt 



(13) 

Since r varies very slowly, we have ignored the second 
order time derivative of r 2 . Now let us take irr 2 = A and 
■nrl = A . Then, 

dA A virdA I 1 4ir(p-p a ) 

dt A rj dt A A mn a r] 

Solving the above equation, assuming that p is indepen- 
dent of r, we get 

A A A ut: vttA 4tt{ P - p a )t 

A y A J A Q i]A Q r)A\ mn a r]A 

(15) 
where c is a constant of integration which must be evalu- 
ated from the boundary condition. Let t = t* be the time 
at which the whole layer is squeezed out. Then for t — t*, 
A=Aq and putting this value in the above equation we 
will have, 

4 + ^)ln(A) + (l-A)(l + ^) = ***L<?-t) 
A r]A A A 7]A mn a r\A a 

(16) 
Now the adiabatic work required to squeeze out one 
monolayer of fluid will be equal to (p — po)A. This 
adiabatic work must be equal to the change in free en- 
ergy due to the process. The total change in free en- 
ergy is (2j si + Vu - 2-f s i + P a)A - V a A, where 7^ 
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FIG. 2: The fit of Eq. (16) to the time sequence of squeezed 
area. The data were obtained from a time sequence picture 
of the squeezed area in the paper of Becker and Mugele. The 
filled squares are the data points and the solid curve is the 
equation fit to the data. 



is the solid-liquid intcrfacial term, Vu is the van der 
Waals interaction term. Equating these two terms we 
get (p-po) = Poa. 

Now Eq. (16) gives the relation between the area of the 
squeezed layer as a function of time. From experiments 
Po, Aq, t* and the three dimensional density, all are 
known. So if we can take a time trace of the area of the 
squeezed monolayer, then by fitting equation (16) to that 
data with 77 and v as two fitting parameters, in principle, 
we can have a value of both 77 and v. So calculating 77 
and v is a matter of having good data. 

Now such good data is lacking in the literature. So, we 
took a time-sequence of pictures of the squeezing out of 
one mono-layer of OMCTS in the paper by Becker and 
Mugele J2Q|. We obtained the area A of the squeezed 
layer as well as the total area Aq of the mica sheets by an 
area-calculating software. Lacking an absolute estimate 
of the areas, we extracted -%- from the pictures, thus 
obtaining the time-evolution of -%-. We took Pq as the 



maximum applied normal pressure Po ~ 2 x 10 N/m 2 
and Aq = 7 x 10 -9 m 2 . These are the typical values of Po 
and A jj, |2{J, UM use(1 in the experiments. 

With the values of the parameters quoted above we 
fitted Eq. (16) to the obtained data, shown in Fig. 2. The 
fit was not very good, since there are two parameters and 
only 8 data-points, but it was as much as could be done 
with the available information. The fit gave the values of 
77 ~ 7 x 10 14 s _1 and v~4x 10 5 cm 2 /s. The value of 77 is 
essentially same as that obtained by Persson and Tosatti. 
For comparison, the kinematic viscosity of OMCTS in 



bulk is 0.02 cm 2 /s. The value of v from the fit is of the 
same order of magnitude as that for confined OMCTS in 
the experiments of 0, [lj| . The experiments of |2J| thus 
indicate an enhancement of viscosity by about ~ 7 orders 
of magnitude. The above result is interesting in at least 
two senses. 

(a) This is an approximate calculation and can not 
be expected to give quantitatively accurate results. In 
this sense it is reassuring that a fit to our theory gives 
a value for v very similar to that in independent exper- 
imental measurements. Crude approximations of layer- 
by-layer flow and p independent of r are made to simplify 
the model, and the dependence of viscosity on the film 
thickness can not be calculated from the model. But at 
molecular thicknesses a liquid film supports a state of 
normal stress and the film thickness adjusts itself to the 
externally applied normal pressure. So, in this sense, the 
thickness and the normal pressure are not two indepen- 
dent parameters. 

(b)From Eq.(16), we see that one of the reasons of high 
increase of the value of v is the strong wall fluid inter- 
action. The viscosity is directly related to the relaxation 
time of a fluid and in the simulation of Scheidler, Kob 
and Binder 7], it is indeed found that the nature of the 
confining walls directly affects the dynamics of the fluid. 

In comparision with the results obtained by Persson 
and Tosatti [8j , we see that the value of 77 is almost same 
as in their calculation. Using the values of the different 
parameters, noted above, r\ ~ 10 14 . But the value of v 
was not obtained by their (incorrect) calculation. The 
model is more interesting than their calculation showed 
it to be, because it tells us that v can be extracted from 
the experiments as well. 

The theoretical approach presented here and, hence, 
the assumption of layer-by-layer flow, can be tested by 
comparison with experiments in which the area squeezed 
out is measured carefully as function of time. In addition, 
different surface treatments should give different degrees 
of anchoring and hence different values of 77, the coeffi- 
cient of surface drag. The trends predicted by Eq. (16) 
with respect to changes in 77 can then be compared with 
experiment. 

The results presented here show that though the model 
is very simple and unable to explain the phenomenon 
completely, it nonetheless demands at least a little at- 
tention. 
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